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Abstract 

We present the results of first-principles calculations on two possible termi- 
nations of the (001) surfaces of SrTi0 3 (STO), BaTi0 3 (BTO), and PbTiOg 
(PTO) perovskite crystals. Atomic structure and the electronic configura- 
tions were calculated for different 2D slabs, both stoichiometric and non- 
stoichiometric, using hybrid (B3PW) exchange-correlation technique and re- 
optimized basis sets of atomic (Gaussian) orbitals. Results are compared with 
previous calculations and available experimental data. The electronic density 
distribution near the surface and covalency effects are discussed in details for 
all three perovskites. 

Key words: SrTi0 3 , BaTi0 3 , PbTi0 3 , single crystal surfaces, surface relaxation, atomic 
and electronic structure, ab initio calculations 



1 Introduction 

The AB0 3 perovskite surfaces are important for many technological applications, 
e.g. STO is widely used as a substrate for epitaxial growth of high- T c superconduct- 
ing thin films, BTO and PTO and their solid solutions are promising for non-volatile 
memory cells, as electro-optical materials, and for piezoelectrical devices [U El H] • 
Due to intensive development and progressive miniaturization of electronic devices, 
the electronic properties and atomic structure of the AB0 3 perovskite thin films 
were extensively studied during the last years. The STO(001) surface structure 
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was analyzed by means of Low Energy Electron Diffraction (LEED) [S], Reflec- 
tion High Energy Electron Diffraction (RHEED), X-ray Photoelectron Spectroscopy 
(XPS) and Ultraviolet Electron Spectroscopy (UPS) |Sj, Medium Energy Ion Scat- 
tering (MEIS) [7 , and Surface X-ray Diffraction (SXRD) [8 . The most recent 
experimental studies on the STO(00f ) include a combination of XPS, LEED, and 
Time-Of-Flight Scattering and Recoiling Spectrometry (TOF-SARS) jHj, and UPS 
and Metastable Impact Electron Spectroscopy (MIES) [TDJ. The BTO and PTO 
surfaces are less studied. The first ab initio calculations of the ABO3 perovskite 
surfaces were presented by Kimura et al. in 1995 [IT] . Since then a number of calcu- 
lations were performed, using different methods, e.g. Linearized Augmented Plane 
Waves (LAPW) [12113, PW ultrasoft-pseudopotential [H H3 HH Ej , Hartree-Fock 
(HF) [THj. Recently, new studies on the STO(OOl) surface relaxation were performed 
using Density Functional Theory (DFT) with plane wave (PW) basis set (DFT-PW 
method) (TH) and Shell Model (SM) simulations with parameters obtained from first- 
principles LAPW method [20 4 .21]. In most of these studies DFT and plane- wave 
basis sets were used (the only exception is HF study [IE])- ^ i s well-known that 
DFT considerably underestimate the band gap. On the other hand, band gap ob- 
tained through the HF calculations is usually an overestimate [22J. To solve this 
problem, the functionals containing "hybrid" of the non-local HF exchange, DFT 
exchange, and Generalized Gradient Approximation (GGA) correlation functionals, 
and known as B3LYP and B3PW (which are quite popular in quantum chemistry 
of molecules) [23 121] could be used. Recently, periodic-structure ab initio hybrid 
calculations were carried out for wide range of crystalline materials [22] , as well as 
for perovskites and their surfaces [2S1 12Z] • In all cases the hybrid functional tech- 
nique shows the best agreement with experimental data for both bulk geometry and 
optical properties of materials under investigation. In this paper we continue our 
own theoretical investigation of the surfaces of perovskite materials. These have 
been carried out using both semi-empirical SM j2H] and ab initio HF and DFT 
methods [2HI EI3 123 EI] • Mostly, these studies are dedicated to STO(OOl) surfaces. 
We studied the effects of different type of Hamiltonians (varied from DFT-LDA to 
HF with posteriori corrections) on surfaces properties and atomic structure. This 
analysis allowed us to choose the B3PW functional as the best for describing the 
bulk properties of all three perovskites under consideration [22] as well as the basic 
surface properties |2"7j . This is why in the present study we use B3PW for compar- 
ing the atomic structure and the electronic properties of the (001) surfaces of three 
similar perovskites, STO, BTO, and PTO. In our simulations, the (001) surfaces of 
perovskites are calculated using a slab model. 

The paper is organized as follows: in Section 2 we present the technical details 
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of calculations and methodology. In Section 3 we discuss the surface structure and 
the electronic properties of the (001) surfaces. Our conclusions are summarized in 
Section 4. 

2 Computational Details 

To perform the first-principles DFT-B3PW calculations, the CRYSTAL '98 computer 
code |22 1321 EH! Oil was use d. This code uses localized Gaussian-type functions 
(GTF) localized at atoms as the basis for an expansion of the crystalline orbitals. 
The ability to calculate the electronic structure of materials within both HF and 
Kohn-Sham (KS) Hamiltonians and implementation of purely 2D slab model with- 
out its artificial repetition along the Z axis, are the main advantages of this code. 
However, in order to employ the LCAO-GTF method, it is desirable to optimize the 
basis sets (BS), which would be suitable for the electronic structure computations. 
Such BS's optimization for all three perovskites is developed and discussed in Ref. 

Unlike the standard basis set [33], we added the polarized O d orbitals, re- 
placed the Ti inner core orbitals by the small-core Hay-Wadt presudopotentials, and 
consistently used the two diffuse s and p Gaussians as the separate basis orbitals on 
the Ti, Ba, Sr and Pb. 

Our calculations were performed using the hybrid exchange-correlation B3PW 
functional involving a hybrid of non-local Fock exact exchange and Becke's three- 
parameter gradient corrected exchange functional jSHj combined with the non-local 
gradient corrected correlation potential by Perdew and Wang E3 EH]- The 
Hay-Wadt small-core ECP's [HHl HOI HI] were adopted for Ti, Sr, and Ba atoms 
[SHI EH E]- The "small-core" ECP's replace only inner core orbitals, but orbitals 
for outer core electrons as well as for valence electrons are calculated self-consistently. 
Light oxygen atoms were left with the full electron BS. The BSs were adopted in 
the following forms: O - 8-411(ld)G (the first shell is of s-type and is a contraction 
of eight Gaussian type functions, then there are three sp-shells and one (i-shell), Ti 
- 411(311d)G, Sr and Ba - 311(ld)G; see Ref. [2E] for more details. 

The reciprocal space integration was performed by the sampling the Brillouin 
zone of the unit cell with the 8x8x1 Pack-Monkhorst net [12], that provides the 
balanced summation in direct and reciprocal lattices 43]. To achieve high accuracy, 
large enough tolerances of 7, 8, 7, 7, 14, (i.e. the calculation of integrals with 
an accuracy of 10~ N ) were chosen for the Coulomb overlap, Coulomb penetration, 
exchange overlap, the first exchange pseudo-overlap, and for the second exchange 
pseudo-overlap respectively [3*2"] . 



3 




a b c 

Figure 1: Schematic view of the modelling for ABO3(001) surfaces: a) AO-terminated 
on both sides, b) Ti02-terminated on both sides, c) asymmetrical termination (AO atop 
and Ti02 on bottom). 

The ABO3(001) surfaces were modelled considering the crystal as a set of crys- 
talline planes perpendicular to the given surface, and cutting out 2D slab of a finite 
thickness, periodic in x-y plane. The symmetrical 7-plane slabs are either AO- or 
TiCVterminated. These are non-stoichiometric with unit cell formulas A4B3O10 
and A3B4O11, respectively. Asymmetrical slab is AO- and Ti02-terminated (from 
each side, respectively), and stoichiometric with a formula A4B4O12 unit cell - are 
shown schematically in Fig. ^ These slabs containing seven planes (symmetric) and 
eight planes (asymmetric) can be considered as thick enough since the convergence 
of calculated slab total energy per ABO3 unit is achieved. This energy differs less 
than 5 • 10~ 4 Hartree for 7- and 9-layered (or 8- and 10-layered for asymmetrical 
termination) slabs for all three perovskites. The symmetric slabs have advantage 
of no dipole moment perpendicular to the slab. However, their non-stoichiometry 
can potentially affect the electronic density distributions and thus atomic displace- 
ments. On the other hand, since the effective atomic charges of Ti and O ions differ 
from the ionic charges +4e, -2e (due to the covalency contribution in Ti-0 chemical 
bonding), the alternating Ti02 and AO planes are slightly charged, which produces 
in the asymmetrical slab certain dipole moment perpendicular to the slab. This 
dipole moment is cancelled by the electronic density redistribution near the surface 
which also can affect the optimized atomic displacements. This is why a critical 
comparison of these two slab models is important for making reliable conclusions. 

In order to compare surface properties of three perovskites, we study here their 
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high-symmetry cubic (Pm3m) phases. The calculated bulk lattice constants (in A) 
are: a = 3.90 for STO, a = 4.01 for BTO, and a = 3.93 for PTO, which demonstrates 
the excellent agreement with experiment • 



3 Results and Discussions 
3.1 Surface Atomic Structure 

In present surface structure simulations we allowed atoms of two outermost surface 
layers to relax along the z axis (by symmetry, surfaces of perfect cubic crystals have 
no forces along the x- and y-axes). Displacements of third layer atoms were found 
negligibly small in our calculations and thus are not treated. The optimization of 
atomic coordinates was done through the slab total energy minimization using our 
own computer code which implements Conjugated Gradients optimization technique 
with a numerical computation of derivatives. 
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Figure 2: Schematic illustration of the relaxation for two outermost surface layers relax- 
ation: a) STO, b) BTO, c) PTO. The view in the [010] direction. Left panels in a-c are 
AO termination, right panels - Ti02- Dashed lines represent positions of the same layers 
in perfect bulk crystals. 



Our calculated atomic displacements are presented in Tableland are schemati- 
cally illustrated in Fig. |21 A comparison with the surface atomic displacements ob- 
tained by other Quantum Mechanical (QM) calculations is also done in Tabled The 
relaxation of surface metal atoms is much larger than that of oxygens which leads 
to a considerable rumpling of the outermost plane. Atoms of the first surface layer 
relax inwards, i.e. towards the bulk. The two exceptions are the top oxygens of STO 
SrO-terminated and PTO Ti02-terminated surfaces (the latter PW-pseudopotential 
calculations are in disagreement however the magnitudes of both displacements 
are relative small, 0.31 and -0.34 per cents of lattice constant, respectively.) The 
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Table 1: Atomic displacements with respect to atomic positions on unrelaxed ABO3(001) 
surfaces (in percent of bulk lattice constant). Symbol A can mean Sr, Ba, or Pb atom. 
Positive displacements are outwards (to the vacuum) , negative displacements mean inward 
the slab center. 
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outward relaxation of all atoms in the second layer is found for all three perovskites 
and both terminations. The displacements obtained for asymmetrically terminated 
slabs are practically the same as for those symmetrically terminated. This means 
that both slab models are reliable for the calculations of the (001) neutral surfaces 
and that 7-8 plane slabs thick enough. 

In order to compare the calculated surface structures with available results ob- 
tained experimentally, the surface rumpling s (the relative displacement of oxygen 
with respect to the metal atom in the surface layer) and the changes in interlayer 
distances Ad±2 and Ac?23 (b 2, and 3 are the numbers of near-surface layers) are 
presented in Table EI Our calculations of the interlayer distances are based on the 
positions of relaxed metal ions (Fig. |2J), which are known to be much stronger elec- 
tron scatters than oxygen ions |5J. The agreement is quite good for all theoretical 
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Table 2: Calculated and experimental surface rumpling s, and relative displacements of 
the three near-surface planes for the AO- and Ti02-terminated surfaces Adij (in percent 
of lattice constant). Results for asymmetrical surfaces are given in brackets. 
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3.9 
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3.76 
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methods, which give the same sign for both the rumpling and change of interlayer 
distances. The amplitude of surface rumpling of SrO-terminated STO is predicted 
much larger than that for TiCVterminated STO surface, whereas the rumpling of 
BTO Ti02-terminated surface is predicted to exceed by a factor of two that for BaO- 
terminated surface. Lastly, PTO demonstrates practically equal rumpling for both 
terminations. From Table El one can see that all surfaces show the reduction of inter- 
layer distance d\ 2 and expansion of (i 2 3- The calculated surface rumpling agrees quite 
well with LEED, RHEED and MEIS experiments jSl El El (which are available for 
STO surfaces only). Theory agrees semi-quantitatively also with the LEED results 
for the A<ii2 and A6?23. However, from Table El is well seen that LEED and RHEED 
experiments contradict each other in the sign of Ag?i2 for SrO-terminated surface 
and A<i 2 3 of Ti02-terminated surface. Another problem is that LEED, RHEED and 
MEIS experiments demonstrate that the topmost oxygen always move outwards the 
surfaces whereas all calculations predict for the Ti02-terminated STO surface that 
oxygen goes inwards. Even more important is the contradiction between the three 
above-mentioned experiments and recent SXRD study [Hj where oxygen atoms are 
predicted to move inwards for both surface terminations reaching rumpling up to 
12.8% for the Ti02 terminated surface! Up to now the reason for such discrepancies 
between the different experimental data is not clear (see [HI EE])- Thus, the disagree- 
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Table 3: Calculated surface energies (in eV per surface cell). Results for previous ab 
initio calculations |lfi| I19 | I2T] are averaged over AO- and Ti02-terminated surfaces. 
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merit between data obtained theoretically and experimentally cannot be analyzed 
until the conflict between different experimental results is resolved. 

The calculated surface energies of the relaxed surfaces, presented in Table were 
computed using the method described in Ref. [2H]- One can see good agreement 
of the surface energies obtained by different methods. The energies calculated for 
AO- and Ti02-terminated surfaces demonstrate a small difference, that means both 
terminations could co-exist. Nevertheless, the energy computed for Ti02-terminated 
STO surface is a little bit larger than that for SrO-termination, in contrast to BTO 
and PTO crystals where Ti02-terminated surface is a little bit energetically more 
favorable. The three (001) surfaces of A 7/ B /y 03 perovskites correspond to "type I" 
stable surfaces revealing the weak-polarity due to partly covalent nature of perovskite 
chemical bonding discussed in the next section. 

3.2 Electronic Charge Redistribution 

We begin discussion of the electronic structure of surfaces with the analysis of charge 
redistribution in near-surface planes. The effective atomic charges (calculated using 
the standard Mulliken population analysis) and dipole atomic moments character- 
izing atomic deformation along the Z axis are presented for all AO-, Ti02- and 
asymmetrically terminated surfaces in Tables EJ El and El respectively. The differ- 
ences in charge densities in the (001) planes in AB0 3 bulk crystals and on the (001) 
surfaces are analyzed in Table 

First of all, note that the effective charges of Sr and Ba are close to the +2e 
formal charges, whereas that of Pb is considerably smaller ( 2e) . Ti and O charges 
are also much smaller than formal charges, similarly to the bulk j2E] which results 
from the Ti-0 covalent bounding. The AO-terminated surfaces of STO and BTO 
show similar behavior. The charges of top layer cations are smaller with respect 
to the relevant bulk charges, the oxygens attract an additional electron charge and 
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Table 4: The calculated Mulliken effective charges and dipole moments for the AO 
termination. Numbers in brackets are deviations from bulk values. Bulk charges (in e); 
STO: Sr = 1.871, Ti = 2.35, O = -1.407, BTO: Ba = 1.795, Ti = 2.364, O = -1.386, PTO: 
Pb = 1.343, Ti = 2.335, O = -1.226 |2SJ- 
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Table 5: For the Ti0 2 
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same as Table 0] 
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Table 6: For the asymmetrical termination. The same as Table 0] 
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0.0067 


-1.414 


-0.0082 


-1.252 


0.0062 






(-0.018) 




(-0.028) 




(-0.026) 




■1 


Ti 


2.352 


0.0026 


2.366 


0.0006 


2.325 


0.0025 






(+0.002) 




(+0.002) 




(-0.010) 






o 2 


-1.408 


-0.0006 


-1.389 


-0.0024 


-1.223 


-0.0323 






(-0.001) 




(-0.003) 




(+0.003) 




5 


A 


1.870 


-0.0042 


1.795 


-0.0074 


1.334 


0.1315 






(-0.001) 




(0.000) 




(-0.009) 









-1.404 


0.0102 


-1.384 


0.0089 


-1.222 


0.0425 






(+0.003) 




(+0.002) 




(+0.004) 




6 


Ti 


2.348 


-0.0089 


2.362 


-0.0084 


2.325 


-0.0115 






(-0.002) 




(-0.002) 




(-0.010) 






o 2 


-1.383 


0.0057 


-1.370 


0.0103 


-1.200 


0.0042 






(+0.024) 




(+0.016) 




(+0.026) 




7 


A 


1.846 


-0.0633 


1.765 


-0.0950 


1.260 


0.0019 






(-0.022) 




(-0.030) 




(-0.083) 









-1.365 


0.0288 


-1.343 


0.0318 


-1.149 


0.0254 






(+0.042) 




(+0.043) 




(+0.077) 




8 


Ti 


2.291 


-0.0748 


2.305 


-0.0823 


2.272 


-0.1508 






(-0.059) 




(-0.059) 




(-0.063) 






o 2 


-1.297 


-0.0493 


-1.278 


-0.0189 


-1.176 


-0.0048 






(+0.110) 




(+0.108) 




(+0.050) 
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Table 7: Calculated charge densities in the (001) planes in the bulk perovskites (in e, 
per Ti02 or AO unit, data are taken from |2H]) and in four top planes of the AO-, Ti02- 
terminated and asymmetrical slabs. Deviations of charge density with respect to the bulk 
are given in brackets. 



Termination 


N 


Unit 


STO 


BTO 


PTO 


Bulk 




AO 


0.464 


0.409 


0.117 






TiOa 


-0.464 


-0.409 


-0.117 


AO 


1 


AO 


0.321 
(-0.143) 


0.278 
(-0.131) 


0.146 
(0.029) 




2 


Ti0 2 


-0.535 
(-0.071) 


-0.457 
(-0.049) 


-0.181 
(-0.064) 




3 


AO 


0.446 
(-0.018) 


0.386 
(-0.023) 


0.096 
(-0.021) 




4 


TiOa 


-0.486 
(-0.022) 


-0.398 
(0.010) 


-0.122 
(-0.005) 


TiOa 


1 


Ti0 2 


-0.334 
(0.130) 


-0.252 
(0.156) 


-0.085 
(0.032) 




2 


AO 


0.490 
(0.026) 


0.422 
(0.013) 


0.104 
(-0.013) 




3 


Ti0 2 


-0.392 
(0.072) 


-0.376 
(0.032) 


-0.078 
(0.039) 




1 


AO 


0.472 
(0.008) 


0.413 
(0.004) 


0.118 
(0.001) 


Asymmetrical 


1 


AO 


0.325 
(-0.139) 


0.281 
(-0.128) 


0.125 
(0.008) 




2 


TiOa 


-0.539 
(-0.075) 


-0.458 
(-0.050) 


-0.166 
(-0.049) 




3 


AO 


0.449 
(-0.015) 


0.386 
(-0.023) 


0.093 
(-0.024) 




4 


TiOa 


-0.464 
(0.000) 


-0.421 
(-0.004) 


-0.121 
(-0.004) 




5 


AO 


0.466 
(0.002) 


0.411 

(0.002) 


0.112 
(-0.005) 




6 


Ti0 2 


-0.418 
(0.046) 


-0.378 
(0.030) 


-0.075 
(0.042) 




7 


AO 


0.481 
(0.017) 


0.422 
(0.013) 


0.111 
(-0.006) 




8 


TiQ 2 


-0.303 
(0.161) 


-0.251 
(0.157) 


-0.080 
(0.037) 
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become more negative. The titanium ions in the second layer demonstrate the slight 
increase of their charges, their oxygens again became more negative due to addition 
electron charge transfer. Changes in atomic charges in deeper layers become very 
small and practically equal zero in the center of the slabs. Unlike STO and BTO, 
the surface oxygen on PbO-terminated surface becomes more positive, Ti in the 
second layer shows practically no changes in the effective charges. 

The charge redistribution in the Ti0 2 -terminated surfaces (Table [5} of all three 
perovskites demonstrates quite similar behavior. All cations in both topmost layers 
demonstrate charge reduction. For surface Ti it is a little bit larger then for A ions 
in subsurface layer. Changes of charges for ions in the asymmetrically terminated 
slabs (Table are practically the same as for the symmetrically terminated AO- 
and Ti0 2 - slab as it should be when slabs are thick enough and the two surfaces do 
not iteract. 

The calculated effective charges of ions in partly-covalent materials usually are 
far from the formal ionic charges due to electron density redistribution caused by 
the covalency effect, what is confirmed by our calculations (see also heading to Ta- 
ble EJ). As a result, the AO and Ti0 2 (001) planes in the bulk perovskites turn 
out to be charged with charge density per unit cell: (Jb{AO) = — crs(Ti0 2 ) (B 
means "Bulk", Table HJ). A half of this charge density comes from Ti0 2 planes to 
each of the two neighboring AO planes. If the formal ionic charges (+4e, -2e, +2e) 
took place, the charge densities of the (001) planes would be zero and we deal with 
the neutral, "type-I surfaces" (according to the generally-accepted classification by 
Tasker jUj). In reality, the charge redistribution makes the ABO 3 (001) surfaces to 
be polar with the dipole moment perpendicular to the surface or "type-Ill polar 
surfaces" , which are potentially unstable. In case of symmetrically terminated slabs 
this dipole moment disappears due to symmetry, but the case of asymmetrical slabs 
is quite of interest. One of possibilities to stabilize the surface is to create surface 
defects |27j . Another option is to add the compensating charge density to the top- 
most layers of surfaces. The charge density changes for all perovskite surfaces with 
respect to the bulk magnitudes are summarized in Table It is well seen, as a result 
of self-consistent calculations, the additional charge densities indeed exist, these are 
localized mainly on two upper layers whereas the central layers practically retain the 
bulk charge density. Such the electron charge density redistribution accompanied 
by atomic displacements on the surfaces allows to compensate the surface dipole 
moment even for the asymmetrically terminated slabs. Indeed, as one can see, the 
sum of changes in charge densities for three topmost layers of the asymmetrically 
terminated slab (on both sides) of all three perovskites approximately equals half 
the charge density of the corresponding bulk plane, which is necessary condition to 
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Table 8: The bond populations for the AO termination (in me, m=milli). Negative bond 
population means atomic repulsion. The corresponding bond populations for the bulk 
perovskites: Ti-0 bond: (STO) 88; (BTO) 100; (PTO) 98; Pb-0 bond: 16. 
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-46 


(V9l 


-38 




-60 




Ti(2) 


78 


Ti(2) 


88 


Ti(2) 


SO 




Sr(3) 


-10 


Ba(3) 


-30 


Pb(3) 


6 




0(3) 


-18 


0(3) 


-34 


0(3) 


-12 


0(3) 


Ti(2) 


86 


Ti(2) 


90 


Ti(2) 


72 




0(3) 


-8 


0(3) 


-6 


0(3) 


-8 




Sr(3) 


-12 


Ba(3) 


-36 


Pb(3) 


24 




Ti(4) 


84 


Ti(4) 


98 


Ti(4) 


96 




0(4) 


-46 


0(4) 


-38 


0(4) 


-54 


0(4) 


Sr(3) 


-10 


Ba(3) 


-31 


Pb(3) 


24 




0(4) 


-8 


0(4) 


-6 


0(4) 


-8 




Ti(4) 


86 


Ti(4) 


98 


Ti(4) 


94 



remove the macroscopic dipole moment These results are in a line with ideas 
of "weak polarity" 0EB]. 

The atomic dipole moments characterize the deformation and polarization along 
the z-axis perpendicular to the surface [32] • These are presented in Tables |H El 
for AO-, Ti02- and asymmetrically terminated surfaces, respectively. On the AO- 
terminated surfaces of all three perovskites the cations have the negative dipole 
moments, directed inwards, to the slab center. The dipole moments of surface Sr, 
Ba and Pb cations on AO-terminated surfaces are surprisingly large, several times 
larger than those of other ions, including the Ti02-terminated surfaces. Oxygens 
of the AO-terminated surfaces of STO and BTO demonstrate small negative dipole 
moments whereas for the PTO the oxygen dipole moment is positive. The polariza- 
tion of second and third layers is pretty small. 

On the Ti02 terminated surfaces polarization of cations has a positive sign, as 
well as for the oxygens in STO and BTO. In contrast, oxygens on PbO surface have 
negative dipole moment. Cations of subsurface layers for all three perovskites reveal 
the positive dipole moments whereas those of oxygens are negative. The asymmetri- 
cally terminated slabs, actually, reproduce the polarization picture obtained for the 
two relevant symmetrically terminated slabs. 

The Mulliken bond populations between atoms in surface layers (which arise 
due to covalency) are presented for the AO-, Ti02- and asymmetrically terminated 
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Table 9: The same as Table |H1 for the TiC>2 termination. 





STO 




BTO 




PTO 




Atom A 


Atom B 




Atom B 




Atom B 




O(l) 


O(l) 


-30 


O(l) 


-24 


O(l) 


-34 




Ti(l) 


111 


Ti(l) 


126 


Ti(l) 


114 




Sr(2) 


-14 


Ba(2) 


-38 


Pb(2) 


42 




0(2) 


-28 


0(2) 


-20 


0(2) 


-42 


0(2) 


Ti(l) 


142 


Ti(l) 


140 


Ti(i) 


162 




0(2) 


2 


0(2) 


2 


0(2) 







Sr(2) 


-8 


Ba(2) 


-30 


Pb(2) 


8 




Ti(3) 


72 


Ti(3) 


90 


Ti(3) 


80 




0(3) 


-.36 


0(3) 


-32 


0(3) 


-36 


0(3) 


Sr(2) 


-4 


Ba(2) 


-24 


Pb(2) 


14 




0(3) 


-42 


0(3) 


-36 


0(3) 


-44 




Ti(3) 


94 


Ti(3) 


106 


Ti(3) 


110 




Sr(4) 


-10 


Ba(4) 


-34 


Pb(4) 


18 




0(4) 


-42 


0(4) 


-36 


0(4) 


-14 


0(4) 


Ti(3) 


92 


Ti(3) 


102 


Ti(3) 


106 




0(4) 


2 


0(4) 


2 


0(4) 


2 




Sr(4) 


-10 


Ba(4) 


-34 


Pb(4) 


14 



slabs in Tables |Hl El and [TUl respectively. The largest effect is observed for PTO 
crystal: the population of the Pb-0 bond on the top layer increases a factor of 
three, as compared to the bulk. The partly covalent nature of Pb-0 bond in lead 
titanate crystal due to hybridization of Pb 6s AO state with the O 2p AO is already 
pronounced in the bulk j2B] , but due to bond shortening (due to surface relaxation) 
its covalency is increased. This effect is also observed on the Ti02 surface, for the 
bond population between surface O and Pb in the second plane. Unlike PTO, in 
STO and BTO there is no indication on the Sr-, Ba- bonding with O atoms. Unlike 
Pb, these cations have effective charges close to formal charge +2e (Table EJ). The 
Ti-0 bonds of all three perovskites on the Ti02 terminated surfaces increase their 
covalency due to bond shortening (caused by surface relaxation) and breaking O 
surface bonds (Table EJ). In line with Ref. [23121], we observe ~ 50% increase 
of the Ti-0 bond covalency on the (001) surface. The asymmetrically terminated 
slabs demonstrate practically the same bond populations as discussed above AO- 
and Ti02-terminated slabs. 

The difference electron density maps calculated with respect to the superposition 
density of spherical A 2+ , Ti 4+ and O 2- ions for STO, BTO and PTO surfaces are 
presented in Fig. El IU and El respectively. These maps demonstrate considerable 
electron charge density redistribution near the perovskite surfaces and are entirely 
consistent with the above-discussed Mulliken charges and bond population analysis. 
For all three perovskites the excess of electron density (the solid iso density curves) is 



14 



Table 10: The same as Table El for the asymmetrical termination. 





STO 




BTO 




PTO 




Atom A 


Atom B 




Atom B 




Atom B 




O(l) 


0(1) 


4 


O(l) 


2 


O(l) 







Sr(l) 


-6 


Ba(l) 


-30 


Pb(l) 


58 




Ti(2) 


70 


Ti(2) 


82 


Ti(2) 


101 




0(2) 


-56 


0(2) 


-60 


0(2) 


-70 


0(2) 


Sr(l) 


-32 


Ba(l) 


-58 


Pb(l) 


48 




0(2) 


-1(5 


0(2) 


-38 


0(2) 


-58 




Ti(2) 


76 


Ti(2) 


88 


Ti(2) 


88 




Sr(3) 


-10 


Ba(3) 


-30 


Pb(3) 


6 




0(3) 


-18 


0(3) 


-32 


0(3) 


-11 


0(3) 


Ti(2) 


86 


Ti(2) 


90 


Ti(2) 


82 




0(3) 


2 


0(3) 


-6 


0(3) 







Sr(3) 


-12 


Ba(3) 


-36 


Pb(3) 


26 




Ti(4) 


86 


Ti(4) 


98 


Ti(4) 


98 




0(4) 


-16 


0(4) 


-38 


0(4) 


-52 


0(4) 


Sr(3) 


-10 


Ba(3) 


-34 


Pb(3) 


14 




0(4) 


-8 


0(4) 


2 


0(4) 


-8 




Ti(4) 


88 


Ti(4) 


98 


Ti(4) 


100 




Sr(5) 


-10 


Ba(5) 


-31 


Pb(5) 


11 




0(5) 


-44 


0(5) 


-38 


0(5) 


-50 


0(5) 


Ti(4) 


88 


Ti(4) 


102 


Ti(4) 


102 




0(5) 


2 


0(5) 


-6 


0(5) 







Sr(5) 


-10 


Ba(5) 


-34 


Pb(5) 


18 




Ti(6) 


88 


Ti(6) 


102 


Ti(6) 


102 




0(6) 


-44 


0(6) 


-36 


0(6) 


-50 


0(6) 


Sr(5) 


-10 


Ba(5) 


-34 


Pb(5) 


24 




0(6) 


-12 


0(6) 


-36 


0(6) 


-44 




Ti(6) 


96 


Ti(6) 


104 


Ti(6) 


110 




Sr(7) 


-4 


Ba(7) 


-24 


Pb(7) 


14 




0(7) 


-36 


0(7) 


-32 


0(7) 


-31 


0(7) 


Ti(6) 


80 


Ti(6) 


90 


Ti(6) 


84 




0(7) 


2 


0(7) 


2 


0(7) 







Sr(7) 


-8 


Ba(7) 


-30 


Pb(7) 


8 




Ti(8) 


132 


Ti(8) 


142 


Ti(8) 


160 




0(8) 


-28 


0(8) 


-20 


0(8) 


-40 


0(8) 


Sr(7) 


-11 


Ba(7) 


-38 


Pb(7) 


42 




Ti(8) 


118 


Ti(8) 
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Ti(8) 
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0(8) 


-32 


0(8) 


-24 


0(8) 


-32 
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STO SrO SID 

' .--I asymmetrical 




Figure 3: The difference electron density maps in the cross section perpendicular to the 
(001) surface ((110) plane) with the AO-, TiC>2 and asymmetrical terminations. Isodensity 
curves are drawn from -0.05 to +0.05 e a.u.~ 3 with an increment of 0.0025 e a.u. -3 . 



16 



BTO BaO BTO 




Figure 4: The BTO difference electron density maps. The same as Fig 
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PTO PbO PTO 




Figure 5: The PTO difference electron density maps. The same as Fig. 
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observed for the Ti-0 bonds, which corresponds to the bond covalency. For all ter- 
minations nearest to the surface Ti-0 bond becomes stronger, but the next nearest 
bond becomes weaker. The A cations on the AO-terminated surfaces demonstrate 
considerable polarization, as it was predicted above from dipole moment calculations 
(Tables EJ El and inj). Nevertheless, the electron density maps demonstrate no clear 
trace of the covalent bonding (zero dot-dashed curves in area between A cations and 
Ti-0 pairs) between A cations and oxygen, even for PTO, despite the Pb-0 bond 
population (Table IBJ) was calculated as 54 me. This means, in reality the covalent 
contribution in Pb-0 bond in PTO is quite weak and plays negligible role in the 
covalency of a whole PTO crystal. Note that the calculated Pb-0 bond population 
considerably depends on the Mulliken population analysis and the Pb atomic BS 
[25] , whose the most diffuse exponent (0.142 bohr -2 ) mainly contributes into the 
increase of covalency. Indeed, recently suggested new approach for the population 
analysis using minimal valence basis of Wannier Type Atomic Orbitals (WTAO) [1H] 
(which are directly connected to the electronic band structure) predicts practicaly 
ionic Pb charge +1.99e for the bulk PTO crystal. Thus, our calculated Pb-0 bond 
population can be interpreted as an increase of the electron attraction between Pb 
and O ions rather a real covalency. 

3.3 Density of States and Band Structures 

The calculated band structures for STO and BTO bulk and surfaces (Fig. |H1 and 
EJ) are quite similar. The band structure for the bulk perovskites was calculated 
using a unit cell which is fourfold extended along the z-axis. Such a supercell is 
similar to the eight-layer "slab" periodically repeated in 3D space and provides also 
the most natural comparison with the surface band structures. In the bulk band 
structure calculations the bands are plotted using the T-X-M-r directions of the 
typical "surface" Brillouin zone (which is simply a square in cubic crystals with the 
standard high symmetry points: T in the center, M the corner and X the center of 
square edge). The upper valence bands (VB) for the STO and BTO bulk are quite 
flat with the top in M point and perfectly flat fragment between M and X points. 
The main contribution into the upper VB comes from O 2p x and 2p y as it is well 
seen from the calculated density of states (DOS) projected onto the corresponding 
atomic orbitals (AOs) (see Fig. 191 and IT2")) . The bottom of lowest conduction band 
(CB) lies at the F point with quite flat fragment between the T and X points and 
consist of Ti 3d threefold degenerated T 2fl level. The optical band gaps for surfaces 
and bulk of all three perovskites as calculated by means of the DFT-B3PW are 
presented in Table HU One can see a good agreement with experiment. We should 
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Table 11: The calculated optical gap (in eV) for the bulk and surface-terminated 
perovskites. The numbers in brackets are from Ref. [2Zj for STO and Ref. [HS] for PTO. 
Both these calculations had no d-orbitals on O, and the later one is done with B3LYP 
functional. The last row contains experimental data. 







STO 






BTO 






PTO 






bulk 


SrO 


TiOa 


asm 


bulk BaO Ti0 2 


asm 


bulk 


PbO 


TiOa 


asm 


Direct gap 


r-r 


3.96 


3.72 


3.95 


3.03 


3.55 3.49 2.96 


2.73 


4.32 


3.58 


3.18 


3.08 




(4.43) 


(4.12) 


(3.78) 








(4.43) 


(3.61) 


(3.77) 




X-X 


4.53 


4.37 


4.04 


4.09 


4.39 4.22 3.63 


3.72 


3.02 


3.79 


3.10 


3.28 




(5.08) 


(4.70) 


(4.38) 








(3.21) 


(3.82) 


(3.12) 




M-M 


5.70 


5.62 


5.17 


4.66 


5.39 5.40 4.17 


4.17 


5.55 


5.37 


5.01 


4.88 




(6.45) 


(5.94) 


(5.04) 








(5.80) 


(6.02) 


(4.89) 




R-R 


6.47 
(7.18) 








6.12 




5.98 








Indirect gap 


x-r 


4.39 


3.55 


3.92 


3.41 


4.20 3.49 3.41 


3.18 


2.87 
(3.18) 


2.96 
(3.03) 


2.98 
(3.12) 


2.78 


M-r 


3.71 


3.30 


3.17 


2.31 


3.60 3.32 2.33 


2.10 


3.66 


3.55 


3.19 


2.96 




(4.23) 


(3.71) 


(3.09) 








(3.85) 


(4.05) 


(2.99) 




R-r 


3.63 
(4.16) 








3.50 




3.66 








LDA-DFT PF (Ref. ITEj) 




1.85 


1.86 


1.13 




1.79 1.80 0.84 




1.54 


1.53 


1.61 




Experiment 






3.75 - direct gap 




3.2 






3.4 







3.25 - indirect gap Ref. 50[ Ref. gl] 

Ref. m 



Table 12: The calculated positions of the valence band top e v and of the conduction band 
bottom e c (in eV) for relaxed and unrelaxed perovskite surface structures. The values in 
brackets are results from Ref. [53] for PTO. The conduction band bottom is in T-point. 
The valence band top lies at the M-point for STO and BTO, but in X-point for PTO. In 
the Ref. [53] the valence band top at the Ti02-terminated surface is in M-point. 



STO BTO PTO 





SrO 


Ti0 2 


BaO Ti0 2 PbO 


TiQ 2 


Unrelaxed structure 


Sc 


1.21 


-2.59 


0.51 -2.86 -1.67 


-2.39 


e v 


-2.48 


-5.22 


-3.01 -5.39 -4.59 


-5.14 


Relaxed structure 


Sc 


-0.50 


-2.78 


0.40 -4.00 -2.11 (-2.13) 


-3.04 (-2.08) 


Sv 


-3.80 


-5.95 


-2.92 -6.33 -5.07 (-5.16) 


-6.02 (5.07) 
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stress here remarkable agreement of the B3PW bulk gap with the experiment for 
STO (3.6 eV vs 3.3 eV). This is in a sharp contrast with the typical HF overestimate 
of the gap and DFT underestimate (e.g. 1.8 eV for STO, BTO [HJEES]). 

The band structure for the SrO-terminated surface demonstrates practically the 
same flat the upper VB as in the bulk STO, with the top of VB at the M point 
and the bottom of CB at the T. The optical band gap for the SrO-terminated 
surface becomes smaller with respect to the band gap of the bulk STO. The narrow 
indirect gap between the T and M points is 3.3 eV, whereas the narrowest gap in 
the bulk is 3.63 eV, i.e. surface gap is reduced by 0.3 eV (see Table [TT] for details). 
Analysis of the Density Of States (DOS) calculated for the SrO-terminated surface 
(Fig. ITU|) demonstrates no contribution of the surface O 2p states into the top of 
VB which mainly consists of 2p AOs of the oxygens from the central plane. The 
main contribution into the CB bottom comes from the Ti 3d which are in the second 
layer. 

The band structure calculated for another, Ti02-terminated surface of STO 
demonstrates less flatness of the VB top, in a comparison with the SrO-termination. 
The indirect optical band gap (M~r) 3.17 eV becomes by fa 0.46 eV narrower as 
compared with the bulk, i.e. twice more reduced as the SrO-terminated surface. 
For the Ti02-terminated STO surface the main contribution into the top of VB is 
made by O 2p x and 2p y AOs which are perpendicular to Ti-O-Ti bridge (see Fig. 
[TTj) . The main contribution to the CB bottom comes from the 3c? AOs of Ti in 
third layer, the energy levels of surface Ti lie a little bit higher in the energy. The 
calculated STO DOS are in a good agreement with MIES and UPS spectra recently 
measured for the Ti02-terminated STO(OOl) surface [HI]- Moreover, our calculated 
position of the VB top for Ti02-terminated STO with respect to the vacuum (5.9 
eV) practically coincides with the experimentally observed value of 5.7 ± 0.2 eV 

(Ref. ma). 

The band structure calculated for the asymmetrically terminated slab demon- 
strates the mixture of band structures obtained for the two symmetrically terminated 
slabs discussed above. The STO band gap becomes narrower (2.31 eV). The VB 
top consists mainly of O 2p from Ti02-terminated slab surface, and the CB bottom 
from 3d AOs of Ti from a subsurface layer (II in Fig. [He)). The split of the upper 
VB (rj 0.8 eV) is well pronounced in asymmetrical STO slab. The band structure 
of the BTO (001) surfaces demonstrates practically the same behavior as STO does 
(Fig. E21 D31 El and Tabled)). Nevertheless, the split of the VB upper band 
is pronounced more for the Ti02-terminated BTO, as compared with the STO sur- 
face. Due to hybridization of Pb 6s and O 2p orbitals in PTO, the calculated band 
structure and DOS slightly differ from those calculated for STO and BTO (Fig. [HI 
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Figure 6: The calculated electronic band structure for STO bulk and surfaces. 
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Figure 7: The calculated electronic band structure for BTO bulk and surfaces. 
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Figure 8: The calculated electronic band structure for PTO bulk and surfaces. 
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Figure 10: Total and projected DOS for the SrO-terminated surface. 
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Figure 11: Total and projected DOS for the STO Ti02-terminated surface. 
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Figure 13: Total and projected DOS for the BaO-terminated surface. 
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Figure 14: Total and projected DOS for the BTO Ti02-terminated surface. 
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EH EE [HI and Table EJ). The narrowest gap of both, bulk and surface band struc- 
tures corresponds to the transition between V and X points of the Brillouin zone. 
In the bulk, the VB top is formed significantly by PB 6s AOs (which also make the 
main contribution to the VB bottom). The bulk CB bottom consist of Ti 3d AOs, 
as in two other perovskites. Surprisingly, the optical band gaps in PbO-terminated 
surfaces are not smaller, as in BTO and STO, but even a little bit increases (up 
to 2.96 eV with respect to 2.87 in the bulk). The VB top in the PbO-terminated 
surface consists of a mixture of Pb 6s and O 2p AOs from a third layer, whereas the 
CB bottom is formed by Ti 3d AOs from a subsurface layer. 

The VB top for the Ti02-terminated PTO(OOl) surface at the X point consists of 
a mixture of the O 2p and Pb 6s AOs from both surface and central layers. Moreover, 
the main contribution comes from the orbitals of central atoms. The CB bottom 
for the Ti02-terminated PTO consists mainly of the Ti 3d AOs from a third layer. 

The band structure calculated for the asymmetrical PTO slab shows a mixture 
of the band structures of the two symmetrical slabs, as it was observed for STO and 
BTO. 

4 Conclusions 

Using hybrid DFT-HF B3PW method, we calculated the surface relaxation and the 
electronic structure of the two possible terminations of the (001) surfaces for STO, 
BTO and PTO perovskite crystals. The data obtained for the surface structure 
are in good agreement with available results of theoretical ab initio calculations 
and with experimental data. The computed relaxed surface energies are quite low, 
1-1.5 eV per unit cell. This indicates that the surfaces with both termination are 
quite stable, in agreement with ideas of a "weak polarity" [3 E3 EE] and existing 
experiments. The analysis of the atomic dipole moments shows that the cations on 
the AO-terminated surfaces reveal the strong electronic polarization along the z-axis. 
The calculated difference electron density maps and bond populations demonstrate 
an increase of the Ti-0 bond covalency near the surfaces, and additionally a weak 
covalency (polarization) of the Pb-0 bond on the PbO terminated surface. 

The observed absence of the surface electronic states in the upper VB for the 
AO-terminated (001) surfaces and considerable ( 0.5 eV) reduction of the gap for the 
Ti02-terminated surface of all three perovskites are important factor for the future 
treatment of the electronic structure of surface defects on perovskite surfaces, as well 
as adsorption and surface diffusion of atoms and small molecules, which is relevant 
for catalysis, fuel cells, and microelectronics. 
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